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In accordance with an old suggestion of Asher Peres (1962), we consider the elec- 
tromagnetic field as fundamental and the metric as a subsidiary field. In following 
up this thought, we formulate Maxwell's theory in a diffeomorphism invariant and 
metric-independent way. The electromagnetic field is then given in terms of the ex- 
citation H = (7i,,T>) and the field strength F = (E,B). Additionally, a local and 
linear "spacetime relation" is assumed between H and F, namely H ~ kF, with 
the constitutive tensor k. The propagation is studied of electromagnetic wave fronts 
(surfaces of discontinuity) with a method of Hadamard. We find a generalized Fres- 
nel equation that is quartic in the wave covector of the wave front. We discuss under 
which conditions the waves propagate along the light cone. Thereby we derive the 
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metric of spacetime, up to a conformal factor, by purely electromagnetic methods. 
PACS numbers: 03.50.De, 04.20.Cv 

Keywords: Classical electrodynamics, premetric axiomatics, differential forms, elec- 
tric/magnetic reciprocity, light cone, metric 



I. INTRODUCTION 

"It is therefore suggested to consider the electromagnetic field as fundamental, and the 
metric field only as a subsidiary quantity." This is a quotation from an article that Asher 
Peres [n| wrote in 1962. Asher's idea was to start with the source-free Maxwell equations 

dH = 0, dF = 0; (1) 

here we have the excitation H = (TC, T>) and the field strength F = (E, B). The equations ((TJ 
are diffeomorphism invariant and independent of the metric gij(x) of spacetime. In other 



words, this is the premetric form [3J, |1 
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22J, |3J| of Maxwell's equations (without 



source). The excitation H and the field strength F are considered to be directly measurable 
quantities, via Maxwellian double place in the case of T> and via the Coulomb force in the 
case of E; analogous procedures exist for H and B, respectively, see 0|. 

Peres then considered the "constitutive" law relating H and F in vacuum, namely 

H = X*F } (2) 

as a definition of the metric; A is a universal constant and * the (metric dependent) Hodge 
star operator. In components the metric is displayed more explicitly, 

= ^V=gg lk g Ji F ke , (3) 

with H l i = ~e l: ' M H^e, where we used the totally antisymmetric Levi-Civita symbol e* jfc£ = 
±1,0. The components of the excitation and field strength 2-forms are given by H = 
Hijdx 1 A dx^ /2 and F = dx % A dx^ /2, respectively. The components H 1 ^ and F^ are 
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assumed to be known and (J3J) has to be resolved with respect to the metric g^. This was 
the program of Peres for deriving the metric from the electromagnetic field. 

Without doubt, Peres is the forerunner of the movement to construct the metric of space- 
time out of purely electromagnetic data. The algebraic method used by Peres in his attempt 
was analyzed in detail by Rubilar j^. Subsequent to Peres Toupin j^Ql and Schonberg 
27 1 proved the existence of a metric in this context. Nowadays we know 4| that the pro- 
gram, as layed out by Peres, can be explicitly implemented. In this paper, which we would 
like to dedicate to Asher Peres on the occasion 70th birthday, we are going to sketch this 
procedure. 

Notation (see ]s]) : We use the formalism of exterior differential forms. We denote the 
frame by e a , with the anholonomic or frame indices a, /?,••• = 0, 1, 2, 3. Decomposed with 
respect to a natural frame <9j, we have e a = e l a <9j, where i,j,--- = 0, 1, 2, 3 are holonomic 
or coordinate indices. The frame e a is the vector basis of the tangent space at each point of 
the 4D spacetime manifold. The symbol J denotes the interior and A the exterior product. 
Both products are metric independent. The coframe $P = efdx^ is dual to the frame, i.e., 



II. PREMETRIC ELECTRODYNAMICS 



Since electric charges are occurring in nature in integer multiples of e/3, here e is the 
elementary charge, they can be counted. Accordingly, if we consider a 3-dimensional (3D) 
volume, we can count how many elementary charges are contained in it. Macroscopically in 
4 dimensions (4D), we can describe the charge density and its flux by the 3- form J that is 
conserved: 

d J = . (4) 

If the global version of charge conservation (0J) is suitably formulated, by de Rham's theorem 
the electric current J turns out to be exact: 

J = dH. (5) 

This is the inhomogeneous Maxwell equation, with the excitation 2-form H. 

In deriving (j3J) and (jSJ), only the ability to define an arbitrary 4D volume is necessary 
together with the counting of electric charges. No distance measurement nor any paral- 
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lei transfer is involved. Therefore it is evident that the inhomogeneous Maxwell equation 
already exists on a 4D manifold without metric and without linear connection. 
With the help of the Lorentz force density 

f a = (e a \F)AJ, (6) 

we can define the electromagnetic field strength 2-form F. If F is integrated over a 2D area, 
it describes the magnetic flux lines piercing through this area. 

In certain situations, inside a superconducting region of a superconductor of type II, 
magnetic flux lines are quantized and can be counted. Here again, the counting of certain 
units, now the magnetic flux quanta, together with the ability to circumscribe a 2D area 
element, is all what is needed to formulate a conservation law. This conservation law of 
magnetic flux we assume to be generally valid: 

dF = . (7) 

Again, no metric and no linear connection of spacetime is involved. As we will discuss below, 
Faraday's induction law is a consequence of (J7J) thereby giving (JJJ) a firm experimental basis. 

Also without using a metric, we can decompose the fields entering the Maxwell equations 
into 1 plus 3 dimensions j^J, 

H = -H Ada + V, J=-jAda+p (8) 

and 

F — E Ada + f\B , A = -ip da + A . (9) 

Here a is the prototype of a time variable. The Lenz factor we put to one in accordance 
with the Lenz rule, f\ = +1, see [^] for details. We also decomposed the potential 1-form A 
that is defined by F = dA. 

Using these decompositions and substituting them into the Maxwell equations © and 
(J7J), we find the conventional form of the Maxwell equations, 

dV = p, \ dB = 0, 

dH = J { dF = I (10) 

T> =dH-j, \ B = -dE, 



cf. Sommerfeld 



291 ] and Scheck 



. The 3D exterior derivative is denoted by d, the time 



derivative by a dot. 
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FIG. 1: Faraday- Schouten pictograms of the electromagnetic field in 3- dimensional space. The 
images of 1-forms are represented by two neighboring planes. The nearer the planes, the stronger 
the 1-form is. The 2-forms are pictured as flux tubes. The thinner the tubes are, the stronger the 
flow is. The difference between a twisted and an untwisted form accounts for the two different types 
of 1- and 2-forms, respectively. Taken from Hehl and Obukhov j^j/. 

We can complete this premetric electrodynamics by substituting (jSj) into ©• Then, by 
some algebra and by some physical arguments, one ends up with an axiom for the energy- 
momentum current 3-form of the electromagnetic field: 



In equations Q to (|TT| . we collected all fundamental relations of classical electrodynam- 
ics. They are all diffeomorphism invariant, that is, completely independent of the coordinates 
used. With the exceptions of (JBJ) and (fTTjl . the equations are also invariant under arbitrary 
linear frame transformations. In contrast, eqs. (JHJ) and (fTTj) . like the frame e a , transform 
linearly, that is, they are covariant under linear frame transformations. 

We started by counting electric charges and magnetic flux lines, used no metric (i.e., 
no gravitational potential) and no linear connection of spacetime, and found a generally 
covariant basis for electrodynamics. In particular, there is no relation whatsoever to the 
Poincare group (also known as inhomogeneous Lorentz group). It simply does not enter at 




(11) 
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that level of our setting up of classical electrodynamics. Classical electrodynamics is not so 
closely related to special relativity as most representations in textbooks want to make us 
believe. 

The present framework can be suitably generalized in order to allow for magnetic charges 
(tantamount of violating magnetic flux conservation, i.e., dF ^ 0), see 0, u|. However, 
since magnetic charges have never been found experimentally, we will not follow this train 
of thought any further. 

In contrast, the violation of electric charge conservation, i.e., dJ ^ 0, would probably 
hit our axiomatic set-up fatally. Firstly, the inhomogeneous Maxwell equation (JHJ) would be 
lost and secondly, the Lorentz force density with a non-conserved charge, is no longer 



expected to qualify as a 
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jona fide defining equation for the electromagnetic field strength 
developed an electrodynamic theory in the framework of 



which one can accommodate a violation of electric charge conservation. They use it as a 
test theory for interpreting corresponding experiments. They propose such experiments in 
order to improve the experimental bounds for charge conservation. 

Since in elementary particle physics all evidence supports electric charge conservation, 
we will stick to this principle and will continue our considerations by assuming the validity 
of (JU) and (J7J), respectively. 

One could ask why these two conservation laws are so "stable" . The absolute dimensions 
of current and field strength are [J] = q and [F] = h/q, respectively, with q ~ dimension 
of charge and h ~ dimension of action. Let us assume that the elementary charge e and 
Planck's constant h are really constants of nature, that is, neither time nor space dependent. 
Then the elementary quanta of charge e/3 and flux nh/e are also constants, and this would 
assure the permanence of charge and flux. The corresponding conservation laws, including 
the premetric Maxwell equations ([10jh are thereby implied. The only problem with this 
argument is that flux occurs in quantized form only under special circumstances. Closely 



related ideas were put forward by Peres 
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Many electromagnetic theories fulfill, indeed, electric charge and magnetic flux conser- 
vation (j3J and (J7|), respectively. Let us quote as examples nonlinear Born-Infeld electro- 
dynamics [lj, the quantum field theoretically motivated pseudo-classical Heisenberg-Euler 
electees ft and , as a more recent case, the axion type electrodynamics of Car- 
roll, Field, and Jackiw Q. Itin Q has shown that the CFJ-electro dynamics can be put 
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in premetric form. We just have to assume the constitutive law H = *F + aF, with the 
metric-dependent Hodge star * and the scalar function a := vt/2, where v is an absolute 
(i.e., prescribed) field and t the time coordinate. A violation of the Lorentz symmetry is 
obtained by postulating this explicit (and noncovariant) time dependence of the constitutive 
tensor. 

As we have seen, premetric electrodynamics turns out to be a useful framework for clas- 
sifying different models of electrodynamics. 



III. LOCAL AND LINEAR SPACETIME RELATION 

The premetric framework, which we discussed so far, is incomplete. As we can read off 
from (|10p. we have 3 + 3 = 6 equations for determining the time evolution of the elec- 
tromagnetic field (H;F). However, the latter is described by the 6 + 6 = 12 components 
(T~C a , Dab = —Dba] E a , B a b = —Bb a ), with a, b, • • • = 1, 2, 3. Thus 6 equations are missing. 

Guided by what we know from ordinary Maxwell-Lorentz electrodynamics, we will as- 
sume that the excitation if is a functional of the field strength F. For vacuum spacetime, 
such a relation is expected to be local, i.e., no integrals are involved linking H(x l ) to F at 
another event (no heredity effects). Furthermore, linearity will be assumed (linear response). 
Accordingly, a local and linear operator k will be postulated to exist, that is, H = k(F). In 
components H = H a /3 , d a A d b /2 etc., this spacetime relation reads, 

H a p = - k Q/ 3 7<5 F y s ■ (12) 

The constitutive tensor of spacetime K a ff & can also be transformed into holonomic coordi- 
nates: K a ff s = e l a ei pCk'ei 5 Kij kl (see Schouten 2^|). Sometimes it is more convenient to 
apply a constitutive tensor with all indices as superscripts. We raise two indices with the 
help of the Levi-Civita symbol and find 

X aM& ■= \e a ^K^ 5 . (13) 

Alternatively, we can introduce a 6-dimensional vector space with the collective indices 
I, K, ■ ■ ■ = 1, . . . , 6 = 01, 02, 03, 23, 31, 12. Then (JT^ can be rewritten as 

Hj = kj K F K = ejMX MK F K ■ (14) 
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All information on the electromagnetic properties of spacetime is encoded in the 36 com- 
ponents of k or Xi respectively. It is straightforward to decompose % irreducibly under the 
linear group GL(4,R). The 6x6 matrix \ 1K decomposes in its symmetric tracefree, its 
antisymmetric, and its trace pieces: 

X IK = (x iIK) ~ ~ e IK XL L ^) + X [m + I e IK XL L , (15) 
36 = 20 © 15 © 1. 

Translated into the 4D formalism, we find 

X HH = (l) x ijkl + e i 3 m[k _ e klm[i pj] + «fci a _ ^ 



■ v ' v 

20, principal 15, skewon 



1, axion 



We indicated the names and the number of independent components in the last equation. 
In conventional Maxwell-Lorentz electrodynamics, only the principal part w^ kl is assumed 
to contribute to the spacetime relation. The tracefree 4x4 matrix ft^ describes the skewon 
part and the axial (or pseudo) scalar a the axion part of \. Again, a transformation from 
anholonomic coordinates can be achieved by the usual rule x^ kl = e% a e> p e fc 7 e l s x a ^ & ■ 
We substitute (fTfij) in the holonomic version of (JTI2J). This yields the spacetime relation 

Htj = l - {1) K tJ kl F kl + 2fi { *F j]k + a Fij . (17) 

Still, the 6D version in 3 + 3 form may be nearer to our intuition. It reads, see Q] for details 
(a,b = 1,2,3), 



V a J 1 Jba V a / \ fib 



(18) 



with the constitutive tensors 



Kj« = ^ Bba , x IK = ^ Va I • (19) 

y A ba v b a J y C a b A ab j 

For material media, such constitutive tensors are well-known, see Lindell et al. Il7| and 

n . 

Mackay and Lakhtakia |18l |. for example. Mackay and Lakhtakia consider a material with 
electric /magnetic cross-terms, i.e., C a b and T> a b are nonvanishing. (However, they require, see 
also (l3j, the axion part to vanish.) It is convenient to write down explicitly the contributions 
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of the principal, the skewon, and the axion parts to the constitutive 3-matrices, cf. (JTfij) : 



A ab = _ e ab _ e abc ^ ^ 

&ab = ^ab + tabc $0°, (21) 

C\ = l\ - $b a - K $ c c ) + a S a b , (22) 

V a b = -fa+tfa'-faW+aft. (23) 



The set of the two symmetric matrices e = e , fi~ b = \x ba plus the traceless matrix ^ a b 
(with 7 C C = 0) comprise the principal part of the constitutive tensor. Usually, e ab is called 
permittivity tensor and fi~ b reciprocal permeability tensor ("impermeability" tensor), since 
they describe the polarizability and magnetizability properties of a medium, respectively 
The cross-term is related to the Fresnel-Fizeau effects. The skewon contributions in (J20J) 
and (|2"Tj) are responsible for the electric and magnetic Faraday effects, respectively, whereas 
skewon terms in ()22j) and (}2"3*|) describe optical activity. 

IV. PROPAGATION OF WAVES: FRESNEL-HADAMARD APPROACH 

Wave propagation is a very important physical phenomenon in classical field theory. In 
electrodynamics, one usually distinguishes between physical optics and geometric optics. We 
will confine our attention to the latter case. The appropriate formalism is then provided by 
the Fresnel-Hadamard approach in which a wave is described in terms of the propagation of 
a discontinuity of the field. Let us define the surface of discontinuity S locally by a function 
$ such that $ = const on S. As usual, we denote by [J 7 ] (x) the discontinuity of a function 
J- across S, and q := is the wave covector. Then for an ordinary wave, the geometric 
Hadamard conditions are satisfied across S: 

[H] = , [dH] = q A h, (24) 
[F] = 0, [dF]=qAf. (25) 

The 2-forms h and / describe the jumps of the derivatives of the electromagnetic fields. 
Using Maxwell's equations (0) and (JTj), we find 

qAh = 0, gA/ = 0. (26) 

The latter equation can evidently be solved by / = q A a, with an arbitrary covector a. Now 
we use the spacetime relation H = k(F). The corresponding relation for the jump 2-forms 
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reads h = n(f) = n(f) + af, where we separated the axion piece and denoted the rest of 
the constitutive relation (|T7|) by the k. Substituting this into we finally obtain 

q A h = q A Tz(q A a) = . (27) 

The last equation can be considered as a kind of Fourier transform of the tensor analytical 
formula dp(j^%As) = , see Post Q, Eq.(9.40). 



A. Generalized Presnel equation and Tamm-Rubilar tensor 

Now we are in a position to derive the Fresnel equation for the wave covector q. As 
a preliminary step we recall that a 3-form in 4D has four components. Consequently, we 
can conveniently recast the 3-form equation (|27|) into an equivalent 4-form equation by 
multiplication with d a = e^dx 1 . Recalling that a = apfi 13 , we find 

q A lT A K(q A ^) a/3 = qA$ a A -^K^q^a^ Af = 0. (28) 

This algebraic system for the covector a obviously admits the gauge freedom ap — > ap + qpip 
with an arbitrary function ip. In order to deal with this problem, we choose the first leg of 
the local anholonomic coframe as i?° = q. This can always be done without restricting the 
generality of our consideration. Then (|28|) reduces to 

iv%i?"Af AfAf = ^ acd Kcd> b a b = (29) 

or, equivalently, to the system of three algebraic equations for the three unknown components 

W ab a b = 6 with W ab :=x 6adb - (30) 

Note that the gauge- dependent czq eventually disappeared. A nontrivial solution of the 
system obtained exists if and only if the corresponding determinant vanishes, 

W := det W ab = y ab ce de fW ad W be W cf = (31) 

or, substituting the components of the 3x3 matrix explicitly, 

U J ± * J ? ~6a0d ~6fe6e ~6c6/ JL n (on\ 

yy - y tabc tdef x x x — u • w z > 
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Quite remarkably, we can rewrite the last equation as a 4D-covariant equation. Indeed, 
because of i a b c = eQafec a ft er some algebra, we find 

w = ^ e mnpq i rstu x mnn x Jpsk x lqtu mnmi ■ (33) 

Here 9 := det(ej a ). Since, by assumption, the 1-forms d a constitute a basis of the cotangent 
space, 9 is always nonvanishing. 

Let us now define 4th-order Tamm-Rubilar tensor density of weight +1, 

<3 ijk \x) ■= ~ e mnpq e rstu X mnr{i X 3[Mk X l)qtu ■ (34) 

In n dimensions a totally symmetric tensor of rank p has ( n+ p _1 ) = ( n ,^lt P ) independent 
components. Accordingly, in four dimensions the TR-tensor (IH4jl has 35 independent com- 
ponents. Because x = X 1 + a e ijfcZ , the total antisymmetry of the Levi-Civita e evidently 
yields G{x) = &(x)- Thus, starting from (JHHj) and discarding the irrelevant nonvanishing 
factor 9 2 , we arrive at the generally covariant 4-dimensional (extended) Fresnel equation 

G m {x)<HWk<li = V- (35) 

This quartic equation in the components g« of the wave covector q was derived from a 
determinant of a 3 x 3 matrix. Apparently, the wave covectors q lie on a quartic Fresnel 
wave surface in general, which, incidentally, is not exactly what we are observing in vacuum 
at the present epoch of our universe. 

B. Decomposing Fresnel equation into time and space pieces 

Physically, the zeroth component component of the wave covector is interpreted as the 
frequency of a "photon", whereas the spatial 3- vector part represents its momentum. Ac- 
cordingly, the physical contents of the Fresnel equation becomes more transparent after 
we carefully separate it into its time and its space pieces. It is convenient to denote the 
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independent components of the TR-tensor (|34|) as follows: 



M 



->0000 



M 



ab 



M 



abc 



M 



I {hell 



<M,A. 



M a := 4 g 000a = -e bcd (A ba A ce C d e + A ab A ec V d ) 
1 



6 £ 00ab = - AW [(C d d ) 2 + (V c c ) 2 - (C c d + V/)(C d c + V, 
+{C d c + V d ){A c{a C b) d + V d {a A b)c ) - C d d A c{a C b) c 
-V c (a A b)c V/ - A dc C {a c V d b) + (A {ab) A dc - A d{a A b)c ) B, 
4 g 0abc = e de(c| [B df {A ab) Vj - V e a A b)f ) 

+B fd (A ab) C{ - A fla C b) e ) + C a f V b) Vj + P/ C b) e C 

r>abcd ef(c \qh\d Y2 



tic 



'0e "e" 9 

Then the Fresnel equation (J3*5)) in decomposed form reads 

g 4 M + g 3 g a M a + q 2 q a q b M ab + g g a g b g c M abc + q a q b q c q d M abcd = . 



(36) 
(37) 



(38) 

(39) 
(40) 



(41) 



V. REDUCING THE QUARTIC WAVE SURFACE TO THE LIGHT CONE 

A. Maxwell-Lorentz electrodynamics in vacuum 

As a first example, let us demonstrate that our formalism yields the correct result for the 
conventional Maxwell-Lorentz spacetime relation with the constitutive tensor density 



X ijkl = 2A 0V /^ i[ V |i] = Aov 7 ^ foV - 9 U 9 Jk ) 



(42) 



Here Ao = \/£o//J>o is the so-called vacuum impedance. A spacetime metric gij is assumed 
on the manifold. Substituting (|4^jl into (jHIJ), we can calculate the corresponding TR-tensor 
density straightforwardly: 



X 3 

Hjkl _ _^3 i /3^ /) (jj„W) _ A 
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-9 {g l3 9 kl + 9 k] 9 U + 9 l3 9 lk ) 



(43) 



As a result, the quartic wave surface (|3~5j) reduces to the usual light cone (twice): 



Q iikl qiWkU = ~ X IV^9 {g ij m)(9 kl (ikqi) = o. 



(44) 



I 



FIG. 2: Clemens Schaefer's drawing 25] of a quartic ray surface for —A = e = e T , /3 = = 1, 
and C = 0, 2? = 0: One ray vector has 4 in general different cuts with the quartic surface. The 
Fresnel equation (|4H|) . which describes the case under consideration, defines a quartic wave surface 
that is dual to the ray surface of this figure. Accordingly, both surfaces encompass the same 
information and look, in fact, fairly similar. Corresponding details are discussed in the Appendix. 

B. Constitutive tensor with vanishing electric/magnetic cross-terms 

The matrices and C b a in the general spacetime relation (jl8j) describe the mixing of 
electric and magnetic fields. Let us consider the case when these cross-terms are absent, 
C ' a = 0, T*b a = 0, whereas A ab and B a b are left asymmetric. Physically, they are interpreted 
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as the generalized permittivity and permeability tensors, respectively. Now the M-tensor 
densities (|36j ) -([40 p simplify appreciably: M a = and M abc = 0, while 

M = detA, (45) 

M ab = ( A (ab) A dc _ A d(a A b)cj ^ ^ ( 4 g) 

M abcd = ±. e ef(* e \9h\b A cd) BgeBhf _ ^ 

The time and space decomposed Fresnel equation (j4*T|) then reduces to 

M ql + M ab q a q b q 2 + M abcd q a q b q c q d = 0. (48) 
This bi-quadratic equation can be solved for the frequency square and yields a Finsler metric, 



see Rubilar 



3, 



, M ab q a q b + VA \ ( 2 , M ab q a q b -VA 



?o + % + 2M 1=°' (49) 



with A := (M ab q a q b ) 2 — 4M M abcd q a q b q c qd- A unique light cone can be recovered, provided 
A = or 

AMM abcd = M {ab M cd)_ ^ 

If this sufficient condition for the existence of a light cone is fulfilled, then the quartic Fresnel 
wave surface reduces to 

/ 2 M ab q a q b \ 2 , ik , 2 , . 

1 go + 2M J ° r \9o V tmk) =0. (51) 

Here the optical metric is introduced with the components g®® t = 2M, g£ b t = M ab . 

If the algebraic condition (J5U|) is not satisfied, the form of the resulting quartic wave 
surface is extremely complicated. For example, even in case when B ab ~ 5 ab and A ab = A ba 
(which describes also the light propagation in an anisotropic purely dielectric crystal), the 
quartic Fresnel surface still looks highly nontrivial, see Fig. 2. 



C. Sufficient condition for unique light cone: electric/magnetic reciprocity 

The algebraic condition (|5()|) imposes a highly nontrivial constraint on the components 
of the constitutive matrices. Although a particular solution (for the homogeneous dielec- 
tric/magnetic medium) can be easily derived, the general solution of (jBTIj) is unknown. There 
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exists, however, another sufficient condition for the reduction of the quartic Fresnel surface 
to a unique light cone which admits a complete solution. Let us put skewon and axion fields 
to zero ft^ = 0, a = 0. Then the spacetime relation (fTTjl contains only the first term on the 
right-hand-side with (fTH|) constructed of the symmetric matrices A ab = A ba , B a b = &ba and 
the traceless matrices T> b = C b a , with C c c = 0. 

Since the times of Maxwell and Heaviside, in the equations of electrodynamics a certain 
symmetry was noticed between the electric and the magnetic quantities and was used in 

n 

theoretical discussions. We formulate electric/magnetic reciprocity as follows |4j: the energy- 
momentum current (jllj) is electric/magnetic reciprocal, i.e., it remains invariant S a — > H a 
under the transformation 

H^(F, F — > — i H , (52) 

with the twisted zero-form (pseudo-scalar function) ( = ((x) of dimension [Q = [H]/[F] = 
q 2 /h = 1/resistance. 

We now require the spacetime relation also to be electric/magnetic reciprocal. Then 

(Fj = Wk/ or - C 2 ^/ = (1) Ki K (1) kk L F l . (53) 

Consequently, the constitutive tensor satisfies a closure relation: 

<»K I K< »K K L = -C 2 6f. (54) 

Therefrom we find ( 2 = — ^Ty(^k 2 ) =: A 2 , and hence (^/t) 2 = —\ 2 1 6 . Thus, symbolically, 
we can write W« = Xy/— 1 6 . With J := \Z-1q, we can then define the duality operator 

*:=J=\^k. (55) 
A 

Besides the closure property J 2 = this operator is also symmetric by construction, 

J ((p) A if) = (f> A J (if)), for all 2-forms (f> and if). As a result, the spacetime relation reads: 

H = X*F with [A] = 1/resistance. (56) 

Using the 3 x 3-matrix parametrization (fT£J) of the constitutive tensor, we can solve the 

C B 

closure relation explicitly. In matrix notation, J = I I . Thus, the closure relation 



T 



A C 

J 2 = — 1 reads: 

C 2 + AB = -1, BC + C T B = 0, CA + AC T = 0. (57) 
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The solution is obtained straightforwardly in terms of A and an arbitrary skew-symmetric 
matrix K = —K T : 



B = -A 



-1 



1 + (KA) 2 , C = KA-\ (58) 



Now we can substitute this solution into (jHf)jl -(|4() )l and find the TR tensor density: 

M = detA, (59) 

M a = Ak a , (60) 

M ab = -2A ab + 6k a % b /detA, (61) 

M abc = 4 (- A (ab k c) + k a k b k c / det A^j /detA, (62) 

M abcd = f A (ab A cd) _ 2 A (ab k c k d) / det A +k a k b k c k d / det A) /detA. (63) 

Here we denote k a := A ab 6b c dK cd . Finally, substituting this into (|41|h we can verify that the 
Fresnel equation reduces to 

(g%q k ) 2 = 0. (64) 



Here the spacetime metric is constructed from the components of the constitutive matrices 
as follows 

9 tk = r-^—r | "V , I (00) 



V- det A y k a - A ab + k a k b / det A , 
It is not difficult to prove that this metric has Lorentzian signature for every matrix A ab . 



VI. DISCUSSION AND CONCLUSION 

The main goal of this paper was to demonstrate that the light cone can be recovered from 
a local and linear spacetime relation of classical electrodynamics. Closure and symmetry of 
this spacetime relation are sufficient conditions that guarantee the reduction of the general 
quartic Fresnel wave surface to a unique light cone. The closure relation alone, without the 
assumption of symmetry, is not sufficient for the recovery of the light cone structure j^j]. 

As a result of the reduction of the quartic Fresnel surface, we find the spacetime met- 
ric with Minkowski (aka Lorentz) signature. This is intimately related to the minus sign 
in the reciprocity transformation (}52|) and the closure relation (|54|) . A plus sign would 
yield the wrong Euclidean signature. Our approach shows that one can treat the du- 
ality operator # as a metricfree predecessor of the Hodge operator * that appears in 
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the standard Maxwell-Lorentz spacetime relation: (duality operator) of Eq. (J55|) — > 
•k (Hodge operator) of Eq. (j2J) . 

Summarizing, the conformal part of the metric, that is, the light cone, naturally emerges 
in our framework from a local and linear spacetime relation. In this sense, the light cone 
(and the spacetime metric) is an electromagnetic construct. 
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VII. APPENDIX: WAVE AND RAY SURFACES 



The drawing in Fig. 2 doesn't depict the wave surface but rather the so-called ray surface 
which is dual to the wave surface. In the particular case of Fig. 2, the propagation in an 
anisotropic dielectric medium is discussed with A ab = —\e ab ) e ab = diag(ei, £2, £3), and 
B a b = c5 ao (with c as speed of light). Then we can immediately verify from (jinj)-(jl7j) that 

M = -6^263/ c 3 , M abcd = c5( ab e cd \ and 

1 



M 



all 



ei(e 2 + e 3 ) 

e 2 (e 1 + e 3 ) • (66) 

\ e 3 (ei + e 2 )j 

As a result, the Fresnel equation for the wave surface (jlSjl can be recast into the simple form 



c 2 (p 



+ 



e 2 q\ 



c 2 (p 



+ 



^4 



c 2 (p 



0. 



(67) 



Here (f = qf + q 2 + q 2 - In crystal optics one usually introduces a ray 4-vector s which is dual 
to the wave covector q, i.e., s\q = 0, see Kiehn et al. [10]. The physical meaning of the ray 
vector is as follows: Its spatial part coincides with the Poynting vector s = E x H, whereas 
its time component describes the ray velocity. Then from ()67|) one can straightforwardly 
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derive the dual equation of the Fresnel ray surface for the components of s = {s°, s 1 , s 2 , s 3 }: 

(^) 2 | (* 2 ) 2 | (* 3 ) 2 _ (68) 

C 2g2 _ ( S 0)2 £x ^ C 2g2 _ ( S 0)2 ^ ^ c 2^ _ ( g 0)2 ^ ^ J 

As we can see, the quartic surfaces (167)1 and ()68j) look pretty much similar, and thus depicting 
one of them in fact gives a good idea about the dual surface. In Fig. 2, we see the ray 
surface with x = s°s 1 /s 2 , y = s°s 2 /s 2 , z = s ^ 3 /^. The vectors & , &, and & represent 
the particular rays related to the corresponding Poynting vectors, whereas D^o, 9t, and W 
represent the wave fronts (propagating along the corresponding wave covectors q) dual to 
them. The Poynting vector s = E x H describes the energy flux density, the wave vector 
q ~ D x B the momentum density. 
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